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ABSTRACT 


Context. Collisional excitation rate coefficients play an important role in the dynamics of energy 
transfer in the interstellar medium. In particular, accurate rotational excitation rates are needed to 
interpret microwave and infrared observations of the interstellar gas for nonlocal thermodynamic 
equilibrium line formation. 

Aims. Theoretical cross sections and rate coefficients for collisional deexcitation of rotationally 
excited HF in the vibrational ground state are reported. 

Methods. The quantum-mechanical close-coupling approach implemented in the nonreactive 
scattering code MOLSCAT was applied in the cross section and rate coefficient calculations on an 
accurate 2D HF-He potential energy surface. Estimates of rate coefficients for H and H 2 colliders 
were obtained from the HF-He collisional data with a reduced-potential scaling approach. 
Results. The calculation of state-to-state rotational quenching cross sections for HF due to He 
with initial rotational levels up to j = 20 were performed for kinetic energies from 10“^ to 
15000 cm“*. State-to-state rate coefficients for temperatures between 0.1 and 3000 K are also 
presented. The comparison of the present results with previous work for lowly-excited rotational 
levels reveals significant differences. In estimating HF-H 2 rate coefficients, the reduced-potential 
method is found to be more reliable than the standard reduced-mass approach. 

Conclusions. The current state-to-state rate coefficient calculations are the most comprehensive 
to date for HF-He collisions. We attribute the differences between previously reported and our 
results to differences in the adopted interaction potential energy surfaces. The new He rate coef¬ 
ficients can be used in a variety of applications. The estimated H 2 and H collision rates can also 
augment the smaller datasets previously developed for H 2 and electrons. 


Key words, molecular processes — molecular data — ISM: molecules 
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1. INTRODUCTION 


Molecular collisions, which are responsible for most of the excitation and reaction processes in¬ 
volving molecules, are important in the interstellar medium (ISM). Collisional excitation and deex¬ 
citation processes compete with radiative transitions in populating molecular levels. In cold envi¬ 
ronments, the important collision partners are H 2 and He because of their large abundances, except 
in photodissociation regions (PDRs) and diffuse gas where collisions with electrons and H can 
become important. Precise laboratory data including collisional deexcitation rate coefficients are 
required for a range of temperatures to interpret the complicated interstellar spectra of molecular 
gas not in local thermodynamic equilibrium (LTE). Because of the complexity and difficulty of 
direct measurements, only limited state-to-state collisional rate coefficients have been measured 
for systems of astrophysical interest (see, for example, [Brunet et al.|2002] l. Therefore, astrophys- 
ical modeling heavily depends on theoretical prediction (e.g., |Flower|2007]|Faure & Lique|2012[ 
[Wiesenfeld & Faure|2013[|Roueff & Lique|2013[|Yang et al.|20T3)|Dubernet et al.|2013| l. 


In this paper, we consider hydrogen fluoride (HF), an interstellar molecule containing a halogen 
element, which was hrst detected in the ISM by Neufeld et al. ( |1997[ ). The reactivity of HF is weak, 
but it may be formed by the exoergic process F- 1 -H 2 —>HF+H. Experimental rate coefficients for this 


reaction were recently reported by Tizniti et al. (2014 1 at temperatures between 10 and 100 K. As a 
result of its structural stability and radiative properties, the HF molecule can be an alternative tracer 
to H 2 in diffuse regions; HF is also the main reservoir of fluorine in the ISM ( [Monje et al.|2011[|van| 
|der Tak |2012a[ ). |Monje et al.| ( |2014| l reported observations of HF in two luminous galaxies, NGC 
253 and NGC 4945, using the Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel 
Space Observatory. With //er.!fc/ief/HIFI, [Neufeld et ^ ( |2010| l detected HF in absorption from the 
ground rovibrational state. Phillips et al.|P010) l presented a detection of the fundamental rotational 
transition of hydrogen fluoride in absorption toward Orion KL using HerschellHWl. The emission 
in the y = 1 ^ 0 rotational transition of HF has been observed in the carbon star envelope IRC 
-1-10216 by |Agundez et al. ( 201 Ij l. Monje et al. ( 2011| l reported the first detection of HF toward a 
high-redshift quasar at z = 2.56, while van der Tak et al.| ( 2012b| l presented observations of the HF 
7=1^0 line in emission towards the Orion Bar. 


The HF-He scattering system has been studied theoretically and experimentally ( |Lovejoy'&] 
[Nesbitt |1990[[Moszynski et al.|1994|[1996[ [Chapman et al.|1997[[Stoecklin et al. [2003[[Reese et[ 

[al.|2005[ [Fajin et aT][2006[ ). The availability of ah initio HF-He potential energy surfaces (PESs) 


has stimulated theoretical studies of HF excitation due to He impact. Lovejoy & Nesbitt (19901 
reported the first experimental study of the near-infrared vibrational-rotational spectra of the HeHF 
and HeDF complexes in a supersonic expansion. The HeHF (HeDF) spectra showed good agree¬ 
ment with the prediction obtained using the Hartree-Fock dispersion (HFD) type rigid-rotor poten¬ 
tial of Rodwell et al. (1981). The spectroscopic data were analyzed and used to probe the isotropic 
and anisotropic intermolecular potentials of the complexes. The two-dimensional (2D) interaction 
potential of [Moszynski et ar| ( [1994| l was developed from ah initio calculations using Symmetry- 
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adapted perturbation theory (SAPT). This SAPT potential is in good agreement with the empirical 
PES of |Lovejoy & Nesbitt | ( |1990| l. All transition frequencies obtained from the bound-state calcu¬ 
lations using the SAPT potential showed excellent agreement with the experimental spectra. The 
SAPT potential has a global minimum for the linear He-HF arrangement and a secondary mini¬ 
mum for the linear He-FH geometry. The accuracy of the SAPT potential was also confirmed by 
agreement between calculated differential and integral cross sections on a slightly modified SAPT 
potential and experimental results ( |Moszynski et al.|1996| l. Another 2D HF-He potential was re¬ 
ported by |Fajin et ah] ( |2006) l using the coupled-cluster method with single and double excitations 
with perturbative triple excitation (CCSD(T)). More recently, a three-dimensional (3D) PES was 
presented by |Stoecklin et al.| ( |2003| l. This PES was computed using the Brueckner coupled-cluster 
method with perturbative triples excitations (BCCD(T)) in the supermolecular approach, and was 
fitted analytically using a kernel Hilbert space interpolation method. This 3D potential was also 
used in close-coupling (CC) calculations of pure rotational excitation of HF in collisions with He 


by Reese et al. (2005 1. Cross sections for transitions for rotational levels up to j - 9 of HF were cal¬ 
culated for collision energies up to 2000 cm“'. Rate coefficients were obtained from 0.1 to 300 K. 
However, the PES of Stoecklin et al. predicts global and local minima with well depths of of 43.70 
and 25.88 cm *, respectively, compared to 39.20 and 35.12 cm ' for the experimentally derived 
PES of |Lovejoy & Nesbitt | ( fT990l l. 

In this work, explicit quantum-mechanical close-coupling scattering calculations of rotational 
quenching of HF in collisions with He to higher levels of rotational excitation were carried out using 
the SAPT potential of |Moszynski et alT] ( |1994| l. The state-to-state rate coefficients are presented for 
a wide range of temperatures (0.1-3000 K), which will aid in modeling rotational spectra of HF 
in various astrophysical and atmospheric environments. The computational method is discussed in 
Sect. 2, and we compare the PESs of Moszynski et al. (1994) and Stoecklin et al. (2003) in Sect. 
3. The results are presented in Sect. 4, while astrophysical applications and scaling approaches for 
estimating HF deexcitation by H 2 and H collisions are described in Sect. 5. 


2. Quantum-scattering calculations 


We adopted the time-independent quantum-mechanical close-coupling (CC) theory developed by 
[Arthurs & Dalgarno~| ( |1963) for the scattering of a linear rigid-rotor by an S -state atom. The state- 
to-state integral cross section for a transition from an initial rotational state j to a final rotational 
state / can be expressed as 


cr 




(Ej) = 




J j J-\- j' 

2(2/+1) 2 Z IV 

/=|y-;l /'=|y-/i 




jj'ir 


-(F/)P 


( 1 ) 


where j and 1 are the rotational angular momentum of the HF molecule and the orbital angular 
momentum of the collision complex, respectively. The total angular momentum J is given by J = 
1 -H j. is an element of the scattering matrix, which is obtained by solving coupled-channel 
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equations and employing the usual boundary conditions, kj = denotes the wave vector 

for the initial channel, Ej is the kinetic energy for the initial channel, and the reduced mass of 
the HF-He system. The total quenching cross section from an initial state j can be obtained by 
summing the state-to-state cross sections <Tj^j>(Ej) over all final / states, where / < j. 

The quantum-scattering code MOLSCAT ( Hutson & Green||1994| l was applied in the close¬ 
coupling calculations. The propagation in R was carried out to 50 A with the coupled-channel 
equations solved using the modified log-derivative Airy propagator ( [Alexander & Manolopoulos] 
|1987[ ). To ensure convergence in the cross-section calculations, at least hve to ten closed channels 
in the basis and a sufficient number of partial waves were included. HF rotational energy levels 


are given in Table 1, which were obtained using the rotational constant Be=20.953 cm ' (Irikura 


20071 and the centrifugal distortion constant £>=0.0021199 cm ' ( Coxon & Hajigeorgiou|[l9^ i. 
The CC calculations were performed for collision energies ranging from 10“^ to 15,000 cm * with 
- 3.3353 u. 

The rate coefficients for rotational transitions can be computed numerically by thermally aver¬ 
aging the corresponding cross sections over a Maxwellian kinetic energy distribution 




(£) = 


_ 8 _ 

TTflfi 


\ 1/2 ^CX 

) 


(Ej) exp(-pEj)dEj, 


( 2 ) 


where T is the temperature, p — (ksT) , and ks is Boltzmann’s constant. 

3. Comparison of potential energy surfaces 

In the rigid-rotor scattering calculations, the interaction potential of HF-He, V(R, 6), was expanded 
in terms of Legendre polynomials. 


V(R,6)^J]^a(R)Pa(cos6), 


( 3 ) 


^=0 


where are Legendre polynomials of order A, R is the distance between the HF center of mass 
and the He atom, and 6 the angle between R and the HF molecular axis. The angular dependence 
of the interaction potential was expanded to Tmax = 20. 

In Fig. [1] the first four components of v^iR), /l=0, 1, 2, and 3 are plotted as a function of R 
for the PESs of Moszynski et al. (1994) and Stoecklin et al. (2003). To obtain a 2D rigid-rotor 
PES, the optimal approach is to average the 3D PES over the ground-state vibrational function 
of the diatomic molecule ( |Eaure et al. [2005] Kalugina et al. 2014| l. However, in this work, we 
mainly used the 2D PES of Moszynski et al. (1994), which was obtained at r^, to calculate the 
rotational quenching cross sections and rate coefficients. Eor the 3D PES of Stoecklin et al., two 
values of the HE bond length, the equilibrium distance r(.=1.7328 ao and the vibrationally averaged 
bond distance for the ground vibrational state ro=1.767 ao ( |Zhang & Zhang |1993[ ), were used to 
compute v^(/?). Comparing the plots, we can see that differences between the PESs of Moszynski 
et al. and Stoecklin et al. are apparent, particularly for components A-1 and 2. In the case of the 
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PES of Stoecklin et al., except for /I = 0, some differences can be seen between the curves for 
R less than 6 ao owing to the different values of the HF bond length. Therefore, one can expect 
discrepancies to arise in scattering calculations performed on the two different PESs. We recall that 


the Moszynski et al. PES agrees very well with the experimentally deduced surface of Lovejoy & 
[Nesbitt | ( fT990l l. 


4. Results and discussion 

4.1. State-to-state and total quenching cross sections 

We performed calculations of the state-to-state quenching cross sections for initial HF rotational 
states of y = 1, 2, ■ ■ ■ , 20 using the PES of Moszynski et al. (1994){^ To evaluate the accuracy 
of the current computed cross section and to compare with the results obtained using the PES of 
Stoecklin et al. (|2003|l, calculations of state-to-state quenching cross sections from j-l and 3 were 


also performed using the PES of Stoecklin et al[^ In the rigid-rotor approximation calculations 
carried out on the PES of Stoecklin et al., the 2D PESs are obtained by fixing the HF bond distance 


at Tg and tq. Correspondingly, rotational constants of HF Be=20.9537 cm ' (Irikura 2007 i and 
Bq - Bg - 1 f2ag-2Q.551Q cm“' were used in the cross-section calculations, where the vibration- 


rotation interaction constant ag = 0.7934 cm“^ (Irikura 2007 i. As examples, the state-to-state 
quenching cross sections from initial states y=l and 3 are presented in Figs.|^and[^ respectively. 
In the case of quenching y = 1 —» 0, Fig. [^illustrates that there are significant differences between 
the cross sections obtained using the PESs of Moszynski et al. and Stoecklin et al. that are due to 
the different stmctures of the two PESs, as shown in Fig.l. For the PES of Stoecklin et al., the cross 
sections calculated using (Bq, tq) agree better with the results obtained using their 3D potential. A 
number of resonances in the cross section, which occur for low collision energies associated with 
the van der Waals wells, demonstrate their sensitivity to the adopted PES. 

The PES of Moszynski et al., used in our calculations, was constructed by fixing the HF bond 
distance at its equilibrium value rg. We used the rotational constant Bg to evaluate the HF rotational 
energy levels (see Table [^. However, to study the effect of the rotational constant on the cross 
sections, we carried out cross-section calculations from initial states j - 1, 3, and 10 using the 
rotational constant Bq. Figures [^ andshow that for the PES of Moszynski et al., the state-to-state 
quenching cross sections from initial j-l and 3 using Bg and Bq are nearly identical. For state-to- 
state quenching cross sections from initial j - 10 (not shown), the differences between the results 
obtained using Bg and Bq are generally lower than 10%. 


* All state-to-state deexcitation cross sections and rate coefficients for HF-He are available on the UGA 
Excitation Database website (www.physast.uga.edu/amdbs/excitation). The rate coefficients are also available 
in the BASECOL ( jPuhernet et al.|2013t and the Leiden Atomic and Molecular Database (LAMDA) ( jSchdierj 
jet al.|2005| l formats. In addition, estimates for HF-H 2 and HF-H rate coefficients obtained by the reduced- 
potential scaling approach, described below, are included in the new LAMDA file along with the original data 
from the LAMDA website. 

^ The rotational deexcitation cross sections displayed in Reese et al. (2005) were found to be discrepant and 
not consistent with rate coefficients given in the same paper as confirmed by T. Stoecklin (priv. communica¬ 
tion). 
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In Fig. 1^ we compare the state-to-state quenching cross section from initial state j - 3. As 
observed in the case of j - 1, the cross sections display resonances in the intermediate energy re¬ 
gion from ~0.01 cm ' to ~10 cm ' due to quasibound levels supported by the attractive part of the 
interaction potential. Importantly, for astrophysical applications, the properties of the resonances 
influence the quenching rate coefficients at low temperatures. In contrast to initial j - the dif¬ 
ference between the cross sections obtained on the two PESs is smaller. In particular, for the PES 
of |Stoecklin et al.| ( |2003| l the results using {Bq, tq) are similar to the results using the 3D potential 
for collision energies higher than 1.0 cm *. As can be seen, the deexcitation process from j - 
is dominated by the Ay = -1 transition. Eurthermore, the computed cross sections show that the 
Ay = / - y = -1 transition dominates the deexcitation for all y, and the cross sections generally 
increase with increasing f , with that for f — 0 being the smallest. 

The total quenching cross section from a given initial level y can be computed by summing over 
all final states f. In Eig.j^the total quenching cross sections from selected initial levels j-2, 4, 6, 
■ ■ ■, 18, and 20 are displayed. Generally, the total quenching cross sections have similar behavior, 
decreasing with y for Ej < 50 cm ', but differences result for small y at high energy that are due 
to a limited number of final exit channels. Each of the cross sections exhibit the behavior predicted 
by the Wigner (1948) threshold law at ultra-low collision energies below ~10“^ cm *, where only 
i-wave scattering contributes and the cross sections vary inversely with the relative velocity. In 
the intermediate energy region, between 0.1 and 10 cm ', the cross sections display scattering 
resonances, but they reveal somewhat different structures depending on the initial rotational state y 
between 0.05 and 1 cm“*. Except for y=2, the total deexcitation cross sections decrease to a global 
minimum near 50 cm"^ 


4.2. State-to-state quenching rate coefficients 


The quenching rate coefficients can be computed by averaging the appropriate cross sections over 
a Maxwell-Boltzmann distribution of collision energy Ej, as given by Eq. The state-to-state 
quenching rate coefficients for initial HE rotational states of y = 1, 2, ■ ■ ■ , 20 were calculated us¬ 
ing Be and the PES of Moszynski et al. (1994). However, to our knowledge, there are no published 
experimental rate coefficients available. Our rate coefficients, computed using the PES of Moszyn¬ 
ski et al., are only compared with the theoretical results of Reese et al. ( |2005 i, which were obtained 
over the limited temperature range of 0.1 to 300 K. As examples, Eigs. m and [^present selected 
state-to-state quenching rate coefficients from initial rotational levels y=l, 3, 5, and 9. Except for 
deexcitation from y = 1, the current results generally follow similar trends with that of Reese et al., 
which were computed using the 3D potential of |Stoecklin et al.| ( |2003[ ). However, at 0.1 K, the rate 
coefficients of Reese et al. are always significantly larger than the current results for all transitions, 
except for y = 1. 

Eor initial state y = 1, Eig.j^a) shows that the current rate coefficients are larger than the results 
of Reese et al. As for initial state y = 3 shown in Eig. |^b), for temperatures above 1 K, the current 
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state-to-state rate coefficients agree reasonably well with that of Reese et al., although our results 
are somewhat larger for the y = 3 —> 1 transition. State-to-state quenching rate coefficients from 
initial states y = 5 and 9 are compared in Figs. |^and|7] respectively, with the results of Reese et 
al. Except for the deexcitation to the final state / = 0 and for 0.1 K, the current rate coefficients 
are smaller than those of Reese et al., similar to what is found for y = 3. For each dominant 
deexcitation transition. Ay = / - y = -1, y=l, 2, • • ■, 9, at a temperature of 50 K, we compare the 
percent differences between our rate coefficients and the results of Reese et al. (2005). In Fig. the 
percent differences are displayed as a function of initial rotational state y. The percent difference is 
near zero for y=3, but the differences vary from 20% to 75% for all othery values. 

For illustration, in Fig. we present the total deexcitation rate coefficients for temperatures 
ranging from 0.1 K to 3000 K for initial states y=10 and 20. Over the whole temperature range 
considered, the rate coefficients generally increase with increasing temperature for all transitions. 
Furthermore, the rate coefficients clearly decrease with increasing |Ay| = |y' - y| with the Ay = -1 
transitions dominant. 


5. Applications 


As highlighted in the Introduction, HF has been observed in both emission and absorption in a 
variety of astronomical environments that may be characterized by diverse physical conditions. As 
a consequence, the rotational levels of HF may be populated by different, or multiple, mechanisms 
leading to spectra differing from FTE. For the UV-irradiated environment of the Orion Bar, |van| 
|der Tak et al.| ( |2012b| l considered electrons and H 2 as possible impactors for inelastic collisional 
excitation. However, given that the Orion Bar is a prototypical photodissociation region (PDR), the 
abundances of both H and He are most likely higher than or comparable to that of H 2 . In many 
PDR enviroments, collisions due to all four colliders may need to be considered. 

For a variety of reasons, it is often the case that excitation rate coefficients for a molecular 
species may only be available for He collisions, as performed in the current work. In such instances, 
a common practice is to estimate H 2 , and occasionally H, collisional excitation rate coefficients us¬ 
ing He data by scaling by the square root of the ratio of the collision systems’ reduced masses, 
here a factor of 1.4 for H 2 . [Walker et al.| ( |2014| l demonstrated on both theoretical and numerical 
grounds that this standard reduced-mass scaling approach is typically invalid. We therfore do not 
recommend here that such an approach be adopted with the present HF-He rate coefficients. Fortu¬ 
nately, explicit HF-H 2 rate coefficients for some rotational quenching transitions of the HF ground 
vibrational state have been computed by Guillon & Stoecklin (2012)l. 


5. 1. Prediction of HF rate coefficients by seating 


Given the availability of HF-H 2 and HF-He collisional data, we can test various scaling methods, in¬ 
cluding the reduced-potential approach introduced in jWalker et al.| ( [2014| l. In the reduced-potential 
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method, the collisional data are scaled by the reduced potentials nx^x according to 




(4) 


where fi is the reduced mass of the HF-Y system, sx is the van der Waals minimum of the HF-Y 
system, and C is a phenomenological exponent. Y = F is usually He, with X - Z typically para-H 2 . 


The reduced-potential and standard reduced-mass approaches are compared in Fig. 10 for the HF 
rotational deexcitation transition y = 1 —> 0. Standard reduced-mass scaling results in an estimate 
for para-H 2 collisions that is a factor of 10 too small, while the reduced-potential method with an 
exponent of C = 1.7 agrees well. 


Given its better performance, the reduced-potential scaling approach of Walker et al. (20141 
was therefore used to predict unknown rate coefficients for para-H 2 and H colliders with HF. We 
adopted the He collision data computed here and the smaller set of para-H 2 rate coefficients for 
j < 5 calculated by |Guillon & Stoecklin| ( p()12| l for T = 10 - 150 K. The rate coefficients are 
very sensitive to the presence of quasibound resonances over this temperature range, which may 
partially be accounted for with the reduced-potential approach since it takes into the account the 
different interaction well-depths. Compared to the reduced-mass scaling technique, where the rate 
coefficients simply scale as the square root of the ratio of the reduced masses, the reduced-potential 
scaling approach offers an improvement for predicted HF rate coefficients, much like the improved 
CO and H 2 O rate coefficient predictions of |Walker et al.| ( |2014| l. While the value of C for CO rate 
coefficients ranged from -0.2 to 1.3 and from 0.5 to 1.2 for H 2 O, it was noted that symmetries 
were involved in obtaining the best value of C. 

The accuracy of the reduced-potential scaling approach can be further tested here with the het- 
eronuclear molecule HF. The normalized root-mean-square deviation (NRMSD), cr„orm^ quantifies 
the residual variance between the calculated H 2 rate coefficients, kcaic, and those scaled from He, 
kscaie, and is given by 


Z (k,cMT)-kcMT))^ 




(5) 


where N is the number of temperature data points and kmax and k^in are the values of the maximum 
and minimum rate coefficients. The resulting NRMSD percentages for the H 2 rate coefficient pre¬ 
dictions for both standard reduced-mass scaling and reduced-potential scaling is given in Table 
and Fig. for the first fifteen transitions of HF at 50 K. When comparing the NRMSD for both 
methods for each transition of HF, the reduced-potential scaling predictions exhibit less residual 
variance in all fifteen transitions with a mean of 14% in its ability to reproduce the explicit HF-H 2 
calculations as computed by |Guillon & Stoecklin| ( |2OT2| l. In many cases, the NRMSD for standard 
reduced-mass scaling exceeds 100%, while the reduced-potential approach gives NRMSD < 35%. 
A linear least-squares analysis was then performed, using the first fifteen transitions of HF, for each 
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/ and the resulting linear functions are plotted in Fig. Except for / = 2, the lines converge 
around |Ay| = 6 and C = -3. For |A= 1, the values of C decrease linearly with / and are listed in 
Table [^for f - 0 - 4. The optimal values for C, valid for all temperatures, were forced to exactly 
reproduce the reduced-potential scaling result for the dominant Aj - -1 transitions. 


To use the reduced-potential approach in astrophysical applications, however, one needs esti¬ 
mates of C when data for the impactor of interest Z, e.g. para-H 2 , are unknown. Using the infor¬ 


mation from Fig. 12 the slope and y-intercept can be obtained for each transition and the value 


of C predicted as given in Fig. 13 and Table 3 for / = 5 - 20. The best prediction for transitions 
with |Ay| > 6 is obtained with C = -3. Rate coefficients for the transitions of H 2 with y < 5 are 
reproduced reasonably well by the reduced-potential approach with C decreasing with increasing 
|Ay| and increasing /. Figures [T0| and [T4| compare rate coefficient estimates using the fit prediction 
of C to direct reduced-potential and reduced-mass values for select transitions. Weak transitions 
are scaled with C - -3, and since these rate coefficients are several orders of magnitude smaller 
than the dominant transitions, larger error in the predicted values is acceptable. 

Given that the available HF-H 2 rate coefficients are limited to y < 5 (|Guillon & Stoecklin 


|2012[ ), we used the reduced-potential scaling method predictions for C and the current HF-He 
data to estimate HF rotational deexcitation for j - 6 - 20 for para-H 2 collisions. Furthermore, as 
there is a complete lack of HF rotational excitation data for H impact, we extended the reduced- 
potential scaling approach to estimate HF-H deexcitation rate coefficients from the current HF-He 


rate coefficients with the trend in C taken from Fig. 13 In using the reduced-potential approach, 
the following parameters were adopted for HF collisions with He, H 2 , and H; /thc = 3.3353 u, 
jUh, = 1.818 u, /in = 0.9596 u, ene = 39.68 cm ' ( Moszynski et al.|[l994 i, = 359.0 cm ' 

( Guillon et al.||20()8 i, and - 100.0 cm ' ( Stark & Werner|fl996 1 . All rate coefficient data are 
available in the FAMDA format, as mentioned in footnote 1. The availability of a complete set 
of HF rotational quenching rate coefficients due to collisions with H 2 , He, H, and e“ will allow 
for detailed modeling of HF rotational emission lines in a variety of environments with a varying 
molecular fraction. However, we caution that the data obtained via reduced-potential scaling are 
approximate, but are reasonable estimates until explicit calculations become available. 


6. Conclusion 


Rate coefficients for the deexcitation of rotational excited HF due to He collisions were computed 
using the close-coupling method and an accurate potential energy surface. The adopted ab initio 
PES agrees well with an experimentally deduced empirical PES. Differences of 75% and larger 
were found with previous HF-He scattering calculations that used a less reliable PES. New rate 
coefficients were obtained for HF rotational levels j - 1-20 for 0.1-3000 K due to He. A recently 
introduced scaling approach was used to estimate rate coefficients for HF-H and missing HF-H 2 
collisional data. 
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Table 1. Computed rotational excitation energies (cm *) of FIR 


j 

Ej 


Ej 

0 

0.00000 

16 

5540.34896 

1 

41.89886 

17 

6210.52398 

2 

125.64498 

18 

6914.70306 

3 

251.13516 

19 

7651.95740 

4 

418.21460 

20 

8421.30660 

5 

626.67690 

21 

9221.71866 

6 

876.26406 

22 

10052.10998 

7 

1166.66648 

23 

10911.34536 

8 

1497.52296 

24 

11798.23800 

9 

1868.42070 

25 

12711.54950 

10 

2278.89530 

26 

13649.98986 

11 

2728.43076 

27 

14612.21748 

12 

3216.45948 

28 

15596.83916 

13 

3742.36226 

29 

16602.41010 

14 

4305.46830 

30 

17627.43390 

15 

4905.05520 




Notes. Ej = Bgjij +V) — Dj^(j + 1)^, where and D are the rotational constant and the centrifugal distortion 
constant of FIF, respectively. 


Table 2. Optimized values of C and their respective normalized root-mean-square deviations (NRMSD) for 
collisional deexcitation transitions of FIF with H 2 scaled via the standard reduced-mass (rm) and reduced- 
potential (rp) methods from FIF-He collisional rate coefficients. 


j / 

C 

NRMSDrm 

NRMSDrp 

1-0 

1.7 

1058.85 

34.90 

2-0 

1.3 

299.75 

22.66 

2-1 

0.6 

161.28 

15.42 

3-0 

-0.2 

53.29 

7.55 

3-1 

-0.1 

30.66 

8.77 

3-2 

-0.2 

62.78 

3.77 

4-0 

-1.6 

69.89 

8.14 

4-1 

-0.3 

42.20 

16.48 

4-2 

-1.3 

73.59 

13.37 

4-3 

-0.8 

90.81 

11.69 

5-0 

-1.2 

54.55 

18.68 

5-1 

-2.0 

61.64 

13.50 

5-2 

-0.5 

44.98 

16.22 

5-3 

-1.3 

66.96 

6.48 

5-4 

-1.2 

86.17 

7.98 
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Table 3. HF reduced-potential fitting parameters. The lower states /, the optimized values of C, the change 
in C (AC), and the slope of the derived linear functions. 


/ 

C 

AC 

Slope 

0 

1.7 

1.1 

-0.94 

1 

0.6 

0.8 

-0.72 

2 

-0.2 

0.6 

-0.56 

3 

-0.8 

0.4 

-0.44 

4 

-1.2 

0.2 

-0.36 

5 

-1.4 

0 

-0.32 

6 

-1.5 

0 

-0.30 

7 

-1.6 

0 

-0.28 

8 

-1.7 

0 

-0.26 

9 

-1.8 

0 

-0.24 

10 

-1.9 

0 

-0.22 

11 

-2 

0 

-0.20 

12 

-2.1 

0 

-0.18 

13 

-2.2 

0 

-0.16 

14 

-2.3 

0 

-0.14 

15 

-2.4 

0 

-0.12 

16 

-2.5 

0 

-0.10 

17 

-2.6 

0 

-0.08 

18 

-2.7 

0 

-0.06 

19 

-2.8 

0 

-0.04 

20 

-2.9 

0 

-0.02 
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R (bohr) 


Fig. 1. Potential expansion terms v^, /1=0, 1, 2, and 3, for the FIF-He PESs of Moszynski et al. (1994) and 
Stoecklin et al. and ro are used to obtain 2D PESs for Stoecklin et al. (2003), while the Moszynski et al. 
PES was constructed for r^. 
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Fig. 2. State-to-state rotational quenching cross sections from initial state y = 1 of HF in collisions with He. 
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Fig. 3. State-to-state rotational quenching cross sections from initial state 7 = 3 of HF in collisions with He. 
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Fig. 4. Total deexcitation cross sections from initial states y=2, 4, 6, 8, 10, 12, 14, 16, 18, and 20 of HF in 
collisions with He with the PES of Moszynski et al. (1994). 
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Fig. 5. State-to-state rotational quenching rate coefficients from initial states j=\ and 3 of FIF that are due 
to collisions with Fie. (a) j = 1, (b) j = 3. Lines: present results with the PES of Moszynski et al. (1994); 
symbols: |Reese et'^ ( |2005] > with the PES of Stoecklin et al. (2003). 
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Fig. 6. State-to-state rotational quenching rate coefficients from initial states j=5 of HF that are due to col¬ 
lisions with He. Lines: present results with the PES of Moszynski et al. (1994); symbols: [Reese et al.| ( |2005) 
with the PES of Stoecklin et al. (2003). 


Article number, page 18 of 26 






















Rate coefficient (cm s 


Benhui Yang et al.: Collisional quenching of highly rotationally excited HF 



Fig. 7. Same as Fig.|^ but for initial rotational state j = 9. 
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Fig. 8. Percent difference of HF-He rate coefficients for the deexcitation of level j for the dominant transition 
Ay = -1 at 50 K between current results using the PES of Moszynski et al. (1994) and the results of Reese et 
al. (2005) with the PES of Stoecklin et al. (2003). 
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Fig. 9. State-to-state rotational deexcitation rate coefficients from initial state j = 10 and j = 20 of FIF in 
collisions with He using the PES of Moszynski et al. (1994). (a) j = 10, (b) j = 20. 
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Fig. 10. HF-para-H 2 I Guillon & Stoecklin||2012 1 and FIF-Fle rate coefficients (current results) compared to 
HF-para-H 2 rate coefficients deduced by various scaling approaches for the y = 1 —> 0 transition. See text for 
discussion. 
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Fig. 11. Normalized root-mean-square deviation (NRMSD) in standard reduced-mass scaling and reduced- 
potential scaling for 15 transitions of HF-para-H 2 , truncated at 100%. 
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Fig. 12. Phenomological constant C as a function of |Ay|. A linear least-squares analysis was performed for 
each / and the resulting linear functions are plotted. Note the convergence of these functions (except for 
/ = 2) around |Ayl = 6 and C = -3. Symbols correspond to C determined for actual HF-para-H 2 and HF-He 
data. 
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Fig. 13. Phenomological constant C as a function of |A j\. As a result of the linear decrease of C with f and 
the convergence around |A j\ = 6 and C = -3, the slope and y-intercept can be obtained for each transition and 
the value of C can be predicted. Symbols and fit lines are for / = 0 - 4, same as in Fig. 12. 
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Fig. 14. Same as Fig. 10, but for the y = 5 —» 1 transition. 
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